%!TEX root = all.tex

\begin{figure}

\subfloat[\bf L$_{2}$ norm]{
  \includegraphics[keepaspectratio=true,scale=0.27]{plots/L2norm_2ndorder_graddis_ice.jpg}
}
\subfloat[\bf L$_{2}$ norm ratio]{
  \includegraphics[keepaspectratio=true,scale=0.27]{plots/L2norm_2ndorder_graddis_ice_ratio.jpg}
}
  
\subfloat[\bf L$_{\infty}$ norm]{
  \includegraphics[keepaspectratio=true,scale=0.27]{plots/Linfnorm_2ndorder_graddis_ice.jpg}
}
\subfloat[\bf L$_{\infty}$  norm ratio]{
  \includegraphics[keepaspectratio=true,scale=0.27]{plots/Linfnorm_2ndorder_graddis_ice_ratio.jpg}  
}

\subfloat[\bf L$_{2}$ and L$_{\infty}$ norms, error ratios and order of accuracy for 3D exponential profile (source terms are from solver's solution)]{

%\begin{table}
\setlength\extrarowheight{3pt}
\newcolumntype{Y}{>{\raggedright\arraybackslash}X}
\begin{tabularx}{6in}{|m{55pt}|m{55pt}|Y|Y|m{55pt}|Y|Y|}
\hline
\bf Grid & \bf L2-norm & \bf Ratio & \bf Observed order of accuracy & \bf Max Error & \bf Ratio & \bf Observed order of accuracy \\
\hline
40x40x40 & 8.3835E-09 & -- & -- & 2.6462E-08 & -- & -- \\
\hline
80x80x80 & 1.8630E-09 & 4.50 & 2.17 & 5.9032E-09 & 4.48 & 2.16 \\
\hline
160x160x160 & 4.5803E-10 & 4.07 & 2.03 & 1.4645E-09 & 4.03 & 2.03 \\
\hline
\end{tabularx}
} %end of subfloat

\caption{L$_{2}$ and L$_{\infty}$ errors and their ratios as functions of time from \ul{\emph{second order accuracy with gradient limiter disabled}} test for exponential 3D profile using source terms from solver's function}


%\end{table}
\end{figure}


%\end{document}